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Abstract 

The scaling behaviour of the edge of the Lee- Yang zeroes in the 
four dimensional Ising model is analyzed. This model is believed to 
belong to the same universality class as the <p\ model which plays a 
central role in relativistic quantum field theory. While in the ther- 
modynamic limit the scaling of the Yang-Lee edge is not modified by 
multiplicative logarithmic corrections, such corrections are manifest 
in the corresponding finite-size formulae. The asymptotic form for 
the density of zeroes which recovers the scaling behaviour of the sus- 
ceptibility and the specific heat in the thermodynamic limit is found 
to exhibit logarithmic corrections too. The density of zeroes for a 
finite-size system is examined both analytically and numerically. 
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1 Zeroes of the Partition Function 



The 4D Ising model is believed to belong to the same universality class as the 
4 model which plays a central role in relativistic quantum field theory. The 
grand canonical partition function for the Ising model (which corresponds to 
the vacuum to vacuum transition amplitude in 4 theory) in the presence of 
an external magnetic field H is 

z N = J_ £ e -«(- j i:«> j (L1) 

N {<M 

where k = (A^T) -1 is the inverse of the Boltzmann constant times the tem- 
perature, J is a coupling constant representing the strength of the intersite 
interaction (set to unity in the following) and N represents the total number 
of sites on the lattice. The value, fa, of the spin at site %, is restricted to 
±1. In its generic form, only nearest neighbour interactions are considered 
and such a link is represented by The sum runs over all M possible 

configurations of the spin field, and the normalization ensures Zn = 1 when 
k — 0. The partition function may be expressed as 

N dN 

Zn= E E p(S,M)e KS e hM , (1.2) 

M=-N S=-dN 

in which d is the dimensionality of the system, h = kH is the reduced external 
magnetic field, and 

N 

S = Y,Mj , M = E0, , (1-3) 

are the configuration energy and magnetization. In S the sum is over the dN 
nearest neighbours or links of a periodic lattice. The spectral density p(S, M) 
denotes the relative weight of configurations having given values of S and M. 
In the absence of an odd external field a second order phase transition occurs 
at a critical value n c of k. The reduced temperature 



is a measure of the distance away from criticality. The partition function Zn 
can be written as an N th degree polynomial in the fugacity z defined by 



e- 2h , (1.5) 



as 

N 

Z N (t,z) = z-^Y.P^) zk > (1-6) 

fc=0 

in which 

dN 

Pk(t)= £ p(S,N-2k)e KS (1.7) 

S=-dN 

is an integrated density. 

That the partition function ( |1.6|) is analytic for finite N establishes that 
no phase transition can occur in a finite-size system. However as N is allowed 
to approach infinity, phase transitions which manifest themselves as points 
of non-analyticity can and do occur. In 1952 Lee and Yang |]| showed 
that the study of the onset of criticality is equivalent to that of the scaling 
behaviour of the zeroes of the partition function. For a finite system, or in 
the thermodynamic limit but in the symmetric phase (t > 0), the zeroes 
in H are strictly complex and the free energy is analytic in a non-vanishing 
neighbourhood of the real axis. As criticality is approached (N — > oo, t — > 0) 
the Lee- Yang zeroes pinch the real H axis, precipitating a phase transition. 
The Lee- Yang theorem states that for the Ising model these zeroes lie on the 
unit circle in the complex fugacity plane (the imaginary axis in the complex 
external field plane). This theorem holds independent of the size, dimension 
and structure of the lattice. 

In the forty years since the ideas of Lee and Yang were presented, there 
has been continual interest in this approach to the problem of phase tran- 
sitions. Analytical progress has included alternative and modified proofs of 
the original circle theorem, extensions of the result to other systems, as well 
as theorems proving that no zeroes can exist in certain regions (see for a 
review) . 

It has been shown rigorously that for isotropic nearest neighbour inter- 
actions, and for t sufficiently positive (the symmetric phase), there exists a 
region around H = which is free from zeroes |3|]. This means there exists 
a gap |Imif| < Hi(t) where the density of zeroes is zero. The free energy is 



3 



analytic in H in the gap and no phase transition can occur (as a function of 
H). The point H = iHi(t), which is a branch point of the partition function, 
is called the Yang-Lee edge || . One expects that this property (the existence 
of a gap) holds in fact for all t > 0. 

Early numerical work on Lee- Yang zeroes involved the exact calculation 
of the density of states p(S, M) and was therefore restricted to very small 
lattice volumes ||. In the 1980's, Monte Carlo histogram approximations to 
the density of states ||, the invention of cluster algorithms |7|, || and multi- 
histogram methods || |10 , ll| all provided boosts to the numerical approach. 
Nonetheless, numerical studies are necessarily limited to finite volume. 

It is, however, the infinite volume limit which is of primary interest. An 
important analytical breakthrough which related such numerical analyses 
in finite volume to the thermodynamic limit came when Itzykson, Pearson 
and Zuber |T2[ connected the concept of partition function zeroes to the 



renormalization group and thereby formulated a finite-size scaling theory 
for these zeroes. Their work applies to dimensions of three or less. This 



was later extended to dimensions of five or more in [DJ. The upper critical 
dimension of the Ising or 4 universality class is d — 4. Above this the 
scaling behaviour of the thermodynamic functions simplifies and the critical 
exponents are exactly those of the mean field theory. Below four dimensions 
the scaling behaviour is of a power-law type. At d = 4 the mean field power- 
law scaling behaviour is modified by multiplicative logarithmic corrections - 
a circumstance intimately related to the expected triviality of the theory flU 



15, 16]. Logarithmic corrections to the finite-size scaling of partition function 
zeroes in four dimensions have recently been identified from a perturbative 
renormalization group analysis backed up by a high precision numerical study 

Recently Salmhofer [18[ has proved the existence of a unique density of 



zeroes in the thermodynamic (N —> oo) limit. The scaling behaviour of the 
Yang-Lee edge in the thermodynamic limit was studied by Abe |T9l and by 



Suzuki 0] in 1967 for Ising models below the upper critical dimension. They 
found asymptotic forms for the density of zeroes and a power-law behaviour 
for the scaling of the edge in the symmetric phase. 

Here we would like to add to this body of knowledge by presenting some 
results on Lee-yang zeroes in the symmetric phase (t > 0) of the four dimen- 
sional Ising model. To this end we find an asymptotic form for the density of 
zeroes in the thermodynamic limit which is sufficient to recover the scaling 
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forms for the specific heat and susceptibility. The finite-size behaviour of 
the location and density of the Lee- Yang zeroes is studied both analytically 
and numerically. 

2 The Density of Lee— Yang Zeroes and the 
Yang— Lee Edge 

According to the Lee- Yang theorem [[[] the zeroes of the partition function 
all lie on the unit circle in the complex fugacity plane. Denoting the (t- 
dependent) position of these zeroes by 

z j (t) = e iB s® , 9je^ , j = l,...,N (2.1) 

the partition function may be written as 

N 

Z N (t, z) = z-*p N {t) II - e^ W ) • (2.2) 

i=l 

The largest coefficient pN{t) plays no role in the following and we henceforth 
set it to unity. The free energy density, 

f N (t,z) = ±\nZ N (t,z) , (2.3) 

can be written as 

1 1 N 

f N (t,z)^--lnz+-YAn(z-e i0 ^) . (2.4) 
The discrete measure dG^ is formally given by 

The t-dependent density of Lee- Yang zeroes on the unit circle in the complex 
z plane is given by and the cumulative density of zeroes Gn is a function 
monotonically increasing in 9 from G(0,t) = to G(2ir,t) = 1. The free 
energy is 

f N (t,z) = --\nz + j g=Q ]n(z-e i6 )dG N (e 7 t) . (2.6) 
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The thermodynamic limit is 

g(9,t) = Jim g N (9,t) , (2.7) 
G(6,t) = lim G N {0,t) , (2.8) 

N—*oo 

f(t,z) = lim f N {t,z) . (2.9) 

The coefficients Pfe(t) of the polynomial (|1.6|) are real and hence g(—6,t) = 
g(8,t). Therefore it is sufficient to consider only the interval < 9 < tc in 
the integrals. The Yang-Lee edge 6 c (t) is defined by 

g(0,t) = for - 6 c (t) < 9 < 6 c (t) . (2.10) 

Integrating ( |2.6| ) by parts gives for the free energy 



/(t,^) = -ln(2cosh(2/ i )+2)- / — — — -G(9,t)d9 . (2.11) 

2 Je c (t) cosh {2h) — cos 9 

The magnetization is then 

<9 f f n sin 9 

= tanh (h) + 2 sinh (2h) / —G(9, t)d9 , (2.12) 

an Je c (t) (cosh (in.) — cos6 l i ) z 

and the zero field susceptibility 

d 2 f\ r n sin9 



x{t) = W)^ 1+i L,w^w Gmde ■ (2 - 13) 

One expects the contribution of small 9 to be dominant Jl9|, |2(| . In particular 
we want to study its contribution singular in t. Expanding the trigonometric 
functions in (|2.13|) (and dropping the constant term), 

X® = 16 r„?&£{l + O{P)}d0 . (2.14) 

J9 c {t) v 

In four dimensions and in the symmetric phase the perturbative renormal- 
ization group gives [21 

X{t)~t- x {-\iit)* . (2.15) 
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A change of variables is introduced via 6 = 6 c x. Then in the critical region 
where t > is sufficiently small 



r^-int) 1 / 3 ~ e c (ty 



2 r /0o(t) G{x6 c ,t) 



dx 



(2.16) 



Following H, H, the upper integral limit can be replaced by infinity near 
criticality. This leads to the requirement that 



t(-lnt)-^ r™G(x6 c ,t) 



e c {ty 



-dx ~ constant. 



x- 



(2.17) 



For fixed t the integral is bounded due to the boundedness of G. The con- 
stancy leads to a differential equation [[II], |2(| for G with the general solution 



G(6,t) =r 1 (-lnt)^ c (t) 2 $ 



e 



Q(x) being an arbitrary function of x with $(| x |< 1) = 0. Then 

6 \ 



--(*). 



(2.18) 



(2.19) 



where $'(x) = d ^}^ ■ 

From ( ^.11| ) (and using the fact that G (0 c ,t) = 0), one gets the specific 
heat 



C v (t) 



d 2 f(t,z) 



dt 2 



h=0 



2 T 9' l ^^-{l + 0(9 2 )}d9 . (2.20) 
Je c (t) dt 2 



Now the cumulative density of zeroes in four dimensions may be found from 
( |2.18D . In four dimensions one expects the power law scaling behaviour char- 
acteristic of dimensions below the upper critical one to be modified by multi- 
plicative logarithmic corrections. Assume therefore that the Yang-Lee edge 
has the scaling behaviour 



6 c (t) = At p (-\nty 
for small t > and with < p < 1. This gives 
d 2 G{6,t) 



(2.21) 



dt 2 



A 2 t 2p ~ 3 (-kit)* 



-2A 



1 + 0| 5T7 



{2(1 - 3p + 2p 2 )$(x) + p{3 - 4p)x$'(x) + p 2 x 2 &'(x)} (2.22) 
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where x = 9/9 c and a prime indicates derivative with respect to x. The 
specific heat is then 



Cv oc t 



2p-3 



-hit): 



-2A 



o 



1 

hit 



57 I(x) 



dx 



x 



where I is some function of x. As t — > (6* c (t) — > 0), one has 



oct 2p " 3 (-lnt)^ 



-2A 



1 + 



1 

hit 



(2.23) 



(2.24) 



in the symmetric phase (t > 0, H — 0) and near criticality. From pertur- 
bation renormalization group analyses it is known |21] that the zero field 
specific heat scales as 

C v (t) ~ (-lnt)s (2.25) 

in four dimensions. Therefore p = | and A = 0. From (|2.21|) the Yang-Lee 
edge in four dimensions scales as 

c (t) ~ • 



(2.26) 



This is the same formula as that yielded by mean field theory |~2], |T^, [20 
The density of Lee- Yang zeroes is given by (|2.19|) as 



g(9,t) (-lnt)3$' 



(2.27) 



in which $' is an unknown function. This form is sufficient to recover the 
singular behaviour of the susceptibility and of the specific heat. 

The behaviour of the zero at 9 = x9 c (for fixed function of t 

(t > 0) is given by (|2.27| ). At fixed t, g(9,t) is an unknown function of 9/9 c . 
Kortman and Griffiths emphasized the study of the density of zeroes close to 
the Yang-Lee edge [p2]] . Using high temperature and high field series they 
concluded that below the upper critical dimension and for a fixed (strictly 
positive) t, the density of zeroes near the edge exhibits a power law behaviour 



g{9,t)~{e-0 e {t)Y 



(2.28) 



In zero dimensions (a single site) a is known to be —1 p3 |. For the 
exactly solvable one dimensional Ising model a = — 1/2 for all t > ffl, i.e., 
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the density of zeroes diverges as the edge is approached. The Ising model 
in the presence of an external field has not been solved in more than one 
dimension. Nonetheless the value of a in two dimensions has been found to 
be — 1/6 by Dhar [[24[] by mapping the two dimensional Ising ferromagnet into 



a solvable model of three dimensional directed animals. Cardy |25] found the 



same result by using the conformal invariance of two dimensional systems at 
the critical point. Using high temperature numerical methods, Kurtze and 
Fisher [|], [2fJ found a = 0.086(15) in three dimensions. It is believed that 
this values hold independent of the lattice parametrization used f§]. For the 
mean field theory a = 1/2 |[22|| . Thus there seems to be a systematic increase 
of a with dimensionality. 

At criticality t = 0, however, the Yang-Lee gap vanishes and one may 
expect the critical exponent o to take on a value different than that in the 
symmetric phase. Now, the density of zeroes is proportional to the disconti- 
nuity in the magnetization M crossing the locus of zeroes [|TJ 

lim M(t, z = re ie ) - lim M(t, z = re ie ) oc g(6, t) . (2.29) 

The infinite volume behaviour of the magnetization below the upper critical 
dimension 

M(t = 0,H) ~ H-s (2.30) 
should be recovered from (|2.28|) at t = and therefore, for d < 4, 

g(9 i t = 0)^9^ . (2.31) 

In four dimensions where 5 = 3, one expects the above formulae to be mod- 
ified by multiplicative logarithmic corrections. There, ( |2.30|) becomes [|2T] 



M(t = 0,H) ~ Hz(-\nH)s . (2.32) 
Therefore, in 4D in the thermodynamic limit 

g (0,t = 0) ~ 05 (-In 0)i . (2.33) 



3 Finite— Size Analysis 

Non-perturbative means of calculating thermodynamic functions in spin 
models are provided by stochastic techniques like Monte Carlo integration. 
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These numerical methods yield exact results subject only to statistical error. 
They are however limited to finite lattices. One has to rely on finite-size scal- 
ing (FSS) extrapolation methods to gain information on the corresponding 
thermodynamic limit. 

Let Piit) represent the value of some thermodynamic quantity P at re- 
duced temperature t on a lattice characterized by a linear extent L. Then, 
if £ is the correlation length, the FSS hypothesis is that Jl7|, [27| 

K$) = f (f^) ' (3,1) 



In four dimensions the scaling behaviour of the correlation length is 28 



Saa(t) ~r*(-lnt)» . (3.2) 

Its scaling with L is [B^] 

&(0)~L(lnL)i . (3.3) 

Therefore in d = 4 the scaling variable should include logarithmic terms [[17]] , 

LQnL)* , ts 

x = -V — i • 3.4 

rl(-int)s 

Let H\ be the position of the Yang-Lee edge in the complex external mag- 
netic field plane in the thermodynamic limit and let Hi(L) be its finite-size 
counterpart (i.e., the position of the lowest lying zero for a system of finite 
linear extent L). From ( |2.26|) , the FSS hypothesis applied to the Yang-Lee 
edge gives 

H X {L) ~ ^ f(x) . (3.5) 

Fixing x (so that when rescaling L, the temperature is also rescaled in such 
a way as to keep x constant), we find 

H X {L) ~ (x- l L- 2 (\nL)-lf f(x) 

~ L- 3 (lnL)~J . (3.6) 

This FSS formula agrees with that derived recently by perturbative renor- 
malization group methods fPT |. 
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The perturbative RG analysis of the finite-size cj)\ model [XT] gives the 
relationship between the magnetization ML(t,H) and external field H at 
reduced temperature t 

H = Cl tM L (\n L) -s + c 2 Ml (In L)" 1 , (3.7) 

where C\ and c 2 are constants. At t — 0, therefore, 

M L (* = 0,#) ~#3(]n£)3 . (3.8) 

For a finite-size system the position of the Yang-Lee edge is not zero at t — 
and the origin of non-vanishing density of zeroes has to be correspondingly 
shifted as in Q2.28 ). One therefore expects the density of zeroes to be 

9L(Hj(L)) ~ (Hj(L) — Hi(L))^ (In L)^ , (3.9) 

where Hj(L) is the position of the j th Lee- Yang zero. Defining the cumulative 
density of zeroes at the j th zero by the fractional total of zeroes up to Hj(L), 

G L {H ] {L)) = 3 -^- , (3.10) 



we find (integrating gi in (|3.9|) to Gl) 
J -I 



Therefore 



(Hj(L) — ifi(L))s (lnL)a . (3.11) 



1 t 

Hj{L) - Ht(L) rsj {-LZ—\ 4 (lnl,)~3 . ( 3 .12) 



L 4 

Eq. ( |3.12| ) gives the j dependence of the lowest lying zeroes as well as recov- 
ering the FSS prediction of (|3.6| ). 

We now compare these FSS results with data obtained for the 4D Ising 
model in a high statistics Monte Carlo calculation. The simulation was done 
with the Swendsen-Wang cluster updating algorithm [^] applied to lattices 
of sizes L 4 for lattice sizes with linear extension L = 8, 12, 16, 20, 24 (details 



of the numerics can be found in 17 



Histogram approximations to the spectral density p(S, M) of ( |1.2| ) were 
determined at zero external field and at various values of k close to the 
pseudocritical value (the value of k at which the zero field specific heat peaks). 
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For the histogram in the magnetization each of the raw histograms in S and 
M were firstly binned in a 256 x 256 array. For each M-bin the corresponding 
S'-subhistograms were then combined to multi-histograms according to [|TO|] . 
In this way, an optimal histogram in M for arbitrary k was obtained. From 
this the partition function may be determined for not too large values of 
(complex) H. 

The critical value of k in four dimensions has been determined to k c = 
0.149703(15) ||17|| . Our data yield only three reliable Lee- Yang zeroes for each 
lattice size. The reason for this is demonstrated in fig.l where the contours 
along which KeZ = and ImZ = (for L = 24 and k = 0.149703) are 
plotted. Because of the magnification of statistical errors far away from the 
simulation point H = these contours fail to cross the imaginary Z axis when 
ImH is large. Thus the zeroes move off the ImH axis and their positions are 
unreliable. The remaining lattices give qualitatively similar pictures. 

Table 1 lists the positions of the first Lee- Yang zeroes (the Lee- Yang 
edge) as obtained from the multihistograms for various k values near k c 
and for all five lattices analyzed. As k increases one expects the zeroes to 
approach the real axis in the thermodynamic limit according to (|2.26| ). Fig. 2 
shows the corresponding behaviour for the finite-size systems considered. At 
n c they should scale according to the FSS formula ( |3.6|) . 

Table 2 lists the positions of the first three Lee- Yang zeroes as obtained 
from the multihistogram at our estimated value for the critical coupling in the 
infinite volume limit, k c = 0.149703. The errors in the quantities calculated 
from the multihistograms were estimated by the jackknife method. The data 
for each lattice size were cut to produce 10 subsamples leading to 10 different 
multihistograms. These 10 different multihistograms give 10 different results 
for the quantities in question, whence the variance and bias were calculated. 

The density of zeroes should behave according to ( |3.9| ) or (equivalently) 
( |3.12|) . The log-log plot of fig.3 gives a slope of 0.778(2). The deviation from 



the exponent 0.75 in ( ^.12[ ) is presumably due to the presence of logarithmic 
corrections. This may be seen in fig.4 where we remove the expected leading 
behaviour: A negative slope is clearly identified. In fact a best fit to all ten 
points gives a slope —0.248(17). The shaded area is bordered by lines of this 
slope. 

We find that both leading power-law scaling behaviour and multiplicative 
logarithmic corrections for the density of zeroes (or equivalently for the dis- 
tance between zeroes) are identified in figs. 3 and 4. This is complementary 
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to our previous analysis in which the scaling behaviour of the actual posi- 
tions of these zeroes was analyzed [[17]]. Both approaches yield quantitative 
agreement with the (perturbative) theoretical predictions. 

4 Conclusions 

The scaling behaviour of the Lee- Yang zeroes and in particular of the Yang- 
Lee edge in four dimensions and in the thermodynamic limit has been exam- 
ined. The asymptotic form for the density of zeroes in the infinite volume 
limit is sufficient to recover the scaling formulae for the specific heat, the 
magnetization and the magnetic susceptibility. This extends the work of 
Abe and Suzuki to the case of four dimensions where mean field power-law 
scaling behaviour is modified by multiplicative logarithmic corrections which 
are linked to the triviality of the theory. An analytical FSS study of the edge 
and the density of zeroes is in good quantitative agreement with a numerical 
analysis in the form of Monte Carlo simulations on finite size lattices. 
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Tables 



Table 1: The positions of the first Lee- Yang zeroes as obtained from the 
multihistograms for all five lattices and near k c . The real part of the zeroes 
is always zero. 







L = 8 


L = 12 


L = 16 


L = 20 


L = 24 




K 


Im#i 


Imtfi 


Im#i 


Im#i 


Im#i 





149600 


0.015281 


0.004511 


0.001958 


0.001047 


0.000637 





149650 


0.015091 


0.004384 


0.001860 


0.000966 


0.000567 





149703 


0.014892 


0.004253 


0.001761 


0.000886 


0.000500 





149750 


0.014718 


0.004140 


0.001677 


0.000820 


0.000447 





149800 


0.014535 


0.004023 


0.001593 


0.000756 


0.000398 





149850 


0.014355 


0.003910 


0.001512 


0.000697 


0.000355 





149900 


0.014177 


0.003800 


0.001437 


0.000644 


0.000319 





149950 


0.014001 


0.003693 


0.001366 


0.000597 


0.000289 





150000 


0.013828 


0.003590 


0.001299 


0.000555 


0.000264 





150050 


0.013657 


0.003491 


0.001237 


0.000518 


0.000244 





150100 


0.013488 


0.003395 


0.001180 


0.000486 


0.000228 





150150 


0.013322 


0.003302 


0.001127 


0.000458 


0.000215 





150200 


0.013159 


0.003213 


0.001077 


0.000433 


0.000204 





150250 


0.012997 


0.003127 


0.001032 


0.000413 


0.000195 





150300 


0.012838 


0.003044 


0.009909 


0.000395 


0.000187 





150350 


0.012682 


0.002965 


0.009530 


0.000379 


0.000180 





150400 


0.012527 


0.002889 


0.009185 


0.000366 


0.000174 



Table 2: The positions of the first three Lee- Yang zeroes as obtained from 
the jackknifed multihistograms at k = 0.149703. The real part of the zeroes 
is always zero. 



L 


Im(#i) 


Im(# 2 ) 


lm(H 3 ) 


8 


0.014892(22) 


0.033057(48) 


0.047357(174) 


12 


0.004253(16) 


0.009426(15) 


0.013349(71) 


16 


0.001761(6) 


0.003905(22) 


0.005388(24) 


20 


0.000886(5) 


0.001970(12) 


0.002743(29) 


24 


0.000500(4) 


0.001106(5) 


0.001541(12) 
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Figures 



Fig 1: Contours along which KeZ = (dotted lines) and ImZ = (full 
lines) (for L = 24 and k = 0.149703). 

Fig 2: The zeroes approach the real axis as k increases; at k c they should 
scale with the lattice size L according to (|3.6|). Here the triangles, circles, 
diamonds, stars and crosses correspond to lattice sizes 8,12,16,20 and 24 re- 
spectively. 

Fig 3: The FSS of the density of zeroes is given by ( p.l2| ). The leading 
power-law behaviour is revealed by a log-log plot. Here the open diamonds 
and triangles correspond to j — 2 and J = 3 respectively. This gives a slope 
0.778(2), the deviation away from 0.75 being due to the presence of logarith- 
mic corrections. 

Fig 4: Data like in fig. 3, but with the leading power-law behaviour removed; 
we clearly identify the negative exponent in the log L behaviour. The shaded 
band indicates the result of a fit giving a slope value of —0.248(17). 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/93 1 1029vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/93 1 1029vl 



